Quenching of High p± Hadron Spectra by Hadronic Interactions 
in Heavy Ion Collisions at RHIC 
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Typically the materialization of high energetic transverse partons to hadronic jets is assumed 
to occur outside the reaction zone in a relativistic heavy ion collision. In contrast, a quantum 
mechanical estimate yields a time on the order of only a few fm/c for building up the hadronic 
wavefunction for jets with typical transverse momenta of p± < 10 GeV as accessible at RHIC 
facilities. The role of possible elastic or inelastic collisions of these high p± particles with the 
bulk of hadrons inside the fireball is addressed by means of an opacity expansion in the number of 
collisions. This analysis shows that the hadronic final state interactions can in principle account for 
the modification of the (moderate) high p± spectrum observed for central collisions at RHIC. 
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I. MOTIVATION AND ESTIMATES 



One of the major goals of the ongoing experi- 
ments at the Relativistic Heavy Ion Collider (RHIC) in 
Brookhaven National Laboratory is to find and probe 
for a temporarily occurring macroscopical state of de- 
confined quark-gluon matter. So far, especially hadronic 
abundancies have been experimentally thoroughly stud- 
ied. Total hadron multiplicity measurements, however, 
do reflect mainly the chemical equilibration and possible 
hadronization processes of the late stages of the system, 
and do not deliver direct information on the early stage 
with maximum energy density and temperature, where a 
deconfined and highly excited state is generally expected 
to be formed. Still, one of the very interesting first re- 
sults of the RHIC experiments is, that one has established 
a significant suppression of moderately high p± hadrons 
produced in central + A collisions compared to rescaled 
peripheral collisions or rescaled (and extrapolated) p + p 
collisions: The results for the lower center-of-mass en- 
ergy of 130 AGeV 0,0 and also the preliminary data for 
200AGeV Q show a suppression factor R(p±) of about 
~ 1/5 for pions with p± sa 5 GeV, and hence state a 
clear hint for (probably various) nuclear medium effects 
at work. 

The most popular explanation for this phenomenon is 
the onset of the occurrence of so called 'jet quenching', 
anticipated in The idea is, that a high energetic par- 
ton moving through a dense coloured (and deconfined) 
medium will loose considerable energy due to collisions 
and induced gluon radiation, and its final fragmentation 
will give rise to particles with considerable lower energies 
Hence, measurements of jets seem to offer a direct 
access to probe the early stage when the deconfined mat- 
ter is very dense. As a present theoretical consensus, 
the dominant mechanism for energy loss is thought to 
be the nonabelian radiation of gluons based on pQCD 
by the energetic partons 0, Q- As a result, the mo- 
menta of the jet partons are attenuated before hadroniz- 



ing. One should internalize the word of caution, that this 
appealing picture of applying perturbative QCD calcula- 
tions is strictly valid only for such high energetic jets 
with E > 10... 20 GeV 0. And indeed, it has been 
found that the suppression pattern of the hard jets (with 
E > 10 GeV) can show some significant dependence on 
the initial properties of the deconfined state §. 

A recent calculation based on the rather involved GLV 
formalism for a finite number of coloured collisions 
does indeed provide a good agreement with the data. 
Still, though, a phenomenological opacity parameter has 
to be adjusted In this study already various other 
effects (the Cronin effect as well as a slight modification 
of the gluon distribution due to shadowing) have been 
phenomenologically incorporated. These effects do partly 
counteract and compete with the pure parton energy loss 
and thus do make a detailed analysis already more deli- 
cate 0, . Contrary, the proposal was made that the 
observed spectra for central collisions show a significant 
m± scaling, being a manifestation of a direct remnant of 
the initial gluons which were liberated from gluon satu- 
rated nuclear distribution functions |l3j . Here any pos- 
sible later interactions of the jets with the surrounding 
medium which might or should alter the distributions are 
completly discarded. 

It has also been raised very recently, that in an ideal 
strong quenching picture the p± spectrum of the high 
momentum hadrons is simply given by early 'surface 
emission' of the outer regions, where roughly half of the 
there produced jets can escape into the vacuum without 
passing through any medium 01 . It has to be verified, 
whether such an idealised scenario of complete absorption 
of energetic partons in matter can actually be achieved. 
The proposal bears resemblance to a similar geometri- 
cal picture for the production of another hard hadronic 
probe, the J/vP, being invoked much earlier for a possi- 
ble explanation of the so called 'anomalous' suppression 
pattern observed in Pb+Pb collisions at CERN-SPS [l5| . 

Typically, it is assumed in all the above refered de- 
scriptions that the partons do exit the collision region as 
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pointlike particles before finally fragmenting into a 'jet' of 
hadrons. We now will point out that this prejudice should 
be of major concern, being invalid for present settings 
at RHIC. The potential magnitude of the hadroniza- 
tion time (or, to be more precise, the time to build up 
the hadronic wavefunction) is based on a relativistic and 
simple quantum mechanical estimate [T(| for light quark 
hadronic states 

t^KER 2 (la) 
or heavy quark meson states 

t^K—R . (lb) 
v m Q 

Here E denotes the energy of the leading 'parton' or jet, 
R states the transversal size of the to be formed hadrons 
and uiq is the mass of a heavy quark. Taking l|la|) with 
the average radius of the pion R^ pa 0.5 fm or (jlb|) substi- 
tuting toq by m p for a p meson and taking R p pa 0.8 fm, 
one has for the formation time a simple understanding 
given as 

t F ps 1 ... 1.2 {Ej GeV) • fm/c . (2) 

Hence, for leading hadrons with moderately high p± < 
10 GeV original point-like jet-partons have established 
already a complete nonperturbative, transversal wave- 
function after traveling a distance in the vacuum of length 
being smaller or equal than 10 fm. Accordingly, the jets 
should, to a large fraction, materialize into hadrons still 
inside the expanding fireball, which has a transversal size 
of roughly 

R{t) w 8fm + 0.5-r, 

where r denotes the local proper time of the longitudi- 
nally (Bjorken like) and transversally expanding system 
with velocity Vt ~ 0.5c. 

It is clear, though, that the parton cannot material- 
ize into hadrons if it still resides in a deconfined phase. 
The dressing with a potential wavefunction and its subse- 
quent fragmentation can only occur outside of deconfined 
matter, i.e. either in the hadronic system or in vacuum. 
The physics of hadronization out of a deconfined plasma 
for fast partons is not understood. Thus the estimate (J2J 
should really be seen as a crude guess. The major fact we 
want to stress is that the final energetic hadrons (with a 
moderately high transverse momentum p±) observed at 
RHIC could well be realized as nearly fully established 
hadrons already inside the late stage hadronic fireball. 
If this is the case, these (pre-)hadrons will interact (by 
collisions) with the bulk hadronic matter making up the 
fireball |l7|. In addition, refering to the phenomenology 
of colour transparency, the colourless objects becoming 
leading hadrons might interact as pre-mesons with the 
surrounding low momentum mesons with a cross section 
of roughly a(t) = <tq ■ {t/tp) as long as their wavefuction 



is not completely being built up [Ti| . This picture of po- 
tential late hadronic final state interaction bears again a 
close similarity to the one applied to explain the J/\V sup- 
pression pattern by either hadronic comovers [l8j or by 
early hadronic string-like excitations as early comovers 

The typical total (elastic plus inelastic) cross section 
of two mesons is given by Co ~ 10 ... 20 mb. Indeed, at 
lower energies for the collisions, there might be strong 
resonance contributions, which give rise to a much larger 
cross section. The density of hadrons in the late fire- 
ball changes from about 1 fm -3 down to 0.1 fm~ 3 . The 
mean free path of the 'fast' hadron is then estimated 
to be A pa 1...10fm. Accordingly, a few collisions 
L/X = 0, 1, 2, 3 . . . should take place. The system is po- 
tentially rather opaque! In return, a strong degradation 
of the momenta of the leading hadrons will result. This 
states the major motivation for the present study: What 
is the potential influence of a specified finite number of 
undergone collisions L/X = 0,1,2,3... of the high mo- 
mentum particles on their spectrum? 

Looking at the available center of mass energy for an 
individual collision for mesons at the same rapidity, even 
for a value of p± = 10 GeV, one gets a y/s < 2 GeV, 
if the target is a pion (assumed being locally) at rest. 
For a p as target, one has an invariant mass of above 
2.5 GeV for particles with a transverse momentum larger 
than 3 GeV. This behaviour is depicted in Fig. ^ I n 
any case, at such lower energies the dominant part of the 
collisions of the mesons or premesons with hadrons of the 
bulk system is nonperturbative and cannot be described 
by pQCD methods. 
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FIG. 1: The CMS energy for three targets at rest {n,K,p, 
depicted by the labels near the curves) and three types of 
projectiles (tt, K, n/p, indicated by the solid, dashed and long 
dashed lines) as function of the momentum of the projectile. 



For the considered transverse momentum region, one 
has either elastic scattering, resonance scattering or also 
inelastic scattering resulting in a few final hadrons. In 
section [n] we explain the description of the collision and 
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introduce the notations of a folding matrix. Since the sur- 
rounding hadronic gas is expected to consist of approxi- 
mately 40% (direct) pions, 10% kaons and 40% rhos, we 
treat in the following a typical inelastic collision within 
the Fritiof scheme |20j with a p-meson as character- 
istic target hadron being at rest. To also stress some 
model dependence, we compare with elastic scattering 
on a 7r as target particle, where the scattering is taken 
simply as isotropic. In section ITTT1 we show various results 
for the transverse momentum spectra at midrapidity for 
various hadronic particles by assuming a definite number 
of undergone collisions. The initial spectra are hereby 
generated by the Pythia code. We end our discussion 
in section IV with a summary and some conclusions. 

II. ENERGY LOSS BY (MULTIPLE) FINAL 
STATE HADRONIC INTERACTIONS 

Matter of interest are the transverse momentum spec- 
tra dNi / dy dp± , where the index i stands for the various 
hadron species and the transverse momentum vector p± 
has to be taken as two dimensional. In all what fol- 
lows, we will assume a Bjorken-like space-time picture 
with azimuthal symmetry, so that the only changes of 
the spectra from the two dimensional transverse momen- 
tum vector are given by its length p± and we therefore 
will skip for reasons of simplicity all dependencies on the 
azimuthal angle. 

In the framework of this assumption we are interested 
in particles lying in the same rapidity window, introduc- 
ing the short hand notation fi(p±) = dNi/dp± for this 
distribution function. (Our notation is different from 
the more common form dNi/p±dp±, which is due to 
practical reasons in the descriptions presented below.) 
In order to separate these distributions of (hard) initial 
hadronic particles of type i from the dominant distribu- 
tions of hadronic particles of type X with moderate or 
low (i.e. 'soft') momenta, say with p± < 2 GeV, we also 
do introduce the distribution functions fx (p± , 4>) ■ Here, 
<f> does not reflect the azimuthal angle in the cylindrical 
shaped Bjorken picture, but denotes the angle between 
the radial transverse direction and its actual momentum 
orientation. These latter distributions do describe the 
motion of the target particles. They are assumed to be 
steeply falling functions in the variable pj_, correspond- 
ing to a rather low temperature in a picture of 'thermal- 
ized' matter, superimposed by some collective, transver- 
sal ('hydrodynamical') flow. 

Within these definitions, one can express the resulting 
spectrum of particles of the kind j after one collision of 
a (hard) particle of kind i with a (soft) target hadron of 
kind X as 

fj(P±) = YlJ d P^ofi(Pu) J dp± x d(f> x (3) 

i 

x ttIx (p± x , <t> x ) gfj [p± ;p± x ,(I> x ;p±) ■ 
l\x 



Here fi(p± a ) nas ^° ^ e interpreted as the initial distri- 
bution of particle kind i, i.e. the distribution one would 
assume without any modifications. Nx denotes the num- 
ber of soft particles of kind X. The "folding matrix" 
9ij (p± :P±x > 4> x expresses the probability to get a 
hadron of species j with momentum p± } if a hadron of 
the species i and the transverse momentum p± collides 
with a target hadron X, which is described by the trans- 
verse momentum component p± x and the orientation <p x . 
Given this general structure with four independent vari- 
ables, gfj is a quite complicated object. It depends on the 
model employed for the collision (see below), and which 
is connected with the target, as indicated by the upper 
index. 

Since the non-transversal component of the target mo- 
mentum, which is encoded in the variable <j) x , also leads 
to scattering in the longitudinal (beam) direction, equa- 
tion |JSj is, strictly speaking, only true if one assumes that 
scattering-out of a specified (tiny) region of rapidity is 
compensated by scattering-in from the next rapidity re- 
gions. This is the situation for the here assumed boost 
invariant Bjorken picture of 1-dimensional longitudinal 
expansion. 

Setting in the following fx(p± x ,<t>x)/Nx = 
8(p± x ) 8((f> x ), we assume for a general estimate 
that the low transverse momentum target is approxi- 
mated by p± ~ 0, which is reasonable, as long as the 
hard momentum scale is clearly separated from the low 
momenta of the bulk hadronic particles. This leads to 
the more simple and self-explanatory folding equation 
for momentum degradation 

fj (PJ.) = 1] / d PU fi (P±o ) 9% (P± ,PX) ■ (4) 

i 

The basic problem is then effectively given by one 
dimensional collisions. Here the reduced folding ma- 
trix (p± n ,p±) has a straightforward interpretation: 
It indicates the probability that out of a given parti- 
cle i with transverse momentum p± one gets a parti- 
cle j with transverse momentum pj_, if the collision is 
taken with a target X at rest, where all this happens 
in the same rapidity interval. E.g. for elastic collisions 
i + X — > i + X (i ^ X), one has the immediate normali- 
sation J °° dp±g(p±_ ,p±) = 1. 

Multiple, subsequent collisions are described by apply- 
ing iteratively eq. f4"|. where in a new iteration step the 
distributions f° have to be replaced by the fi resulting 
from the former step. 

We do remark some words of caution: 
(1) The separation of the transversal distribution into a 
"hard" contribution fi and a "soft" part fx and setting 
the latter to a delta distribution implies some drawback, 
as this separation is artificial at low p± (in the detector 
one would see fi + fi) and thus leads to a "particle pro- 
duction" effect at low values of p±. This is connected 
with the normalization of the folding matrix gfj , what 
one can see easily, if one for example looks at the simple 
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situation of elastic scattering ir + on a 71", i.e. i = ir + and 
X = 7r = +7t° + 7r~), leading to two pions with non 
vanishing transverse momenta. We remark here that one 
has to clearly stay in the discussion above a certain mod- 
erately large transversal momentum scale for an assumed 
clear separation of hard and low momentum scales. 
(2) All our considerations assume, that a fixed num- 
ber of collisions takes place, disregarding the fact, that 
the cross section could vary dramatically with the trans- 
verse momentum. For special collision partners, for 
example p + it, it is known that the cross sec- 

tion a tot {p + [ 7r+ + /2) changes from « 30 mb for 
Pproton = 3 ... 6 GeV to nearly five times this value at 
Pproton ~ 1.8 GeV, if the pion is considerd to be at 
rest. Similarly, the elastic reaction tt + tt — ► p — ► n + n 
also exceeds 100 mb for the projectile momentum p v w 
2 ... 2.5 GeV. This implies that, while assuming that one 
collision appears, at an other value of the transverse mo- 
mentum (in these special cases) five collisons are to be 
awaited. We will comment on this point a little further 
in the conclusions. 



A. Initial distribution 

For the initial distribution ff of individual hadrons i 
at moderately high p± when entering the final hadronic 
stage we do not have direct information. All various 
ideas, as mentioned in the introduction, might contribute 
and already steepen the spectrum. In order to see clearly 
the potential effect of the late hadronic interactions, we 
use individual distributions generated by Pythia v6.2 
22] gauged to individual pp-collisions. 

It is essential to carefully adjust the parameters for the 
model calculations to data provided by the UA1 collabo- 
ration for pp collisions in the energy range s/s > 200 GeV 
and \y\ < 2.5 |23j. Fitting the calculations to the data, 
one gets a best description of the experimental results by 
just tuning one parameter, the intrinsic transverse mo- 
mentum distribution of the partons inside the nucleon, 
described by its averaged value (k±) (yfs). The resulting 
spectra are displayed in Fig. 

In a second step one has to compare the calculated 
spectra with the parametrization (1 + p±/po)~ n used to 
describe the experimental data |23| and extract from this 
fitted parmetrization (po(y/s) and n(y/s)) the character- 
istic value (pj_) With this knowledge one can ex- 
trapolate the calculations of pp collisions down to the 
energy y/s — 130 GeV by adjusting the intrinsic parame- 
ter (k±) to finally meet the experimentally extrapolated 
value (p±)(y/s = 130 GeV) = 0.385 GeV Please 
note that this value is completely disregarded by the as- 
sumed pointwise interpolation of fits to data from differ- 
ent values of v/s for given transverse momentum px, as 
used eg. in |lll24|. 

After having adjusted the microscopic Monte-Carlo 
calculations for pp collisions at y/s = 130 and 200 GeV 
and \y\ < 2.5, theoretical predictions for Au-Au col- 
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FIG. 2: Charged hadron distributions for pp collisions cal- 
culated with Pythia at y/s = 200, 500 GeV compared with 
experimental data [2^|. 
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FIG. 3: Hadron distributions at midrapidity for Au-Au col- 
lisions at y/s = 130 GeV (left) and y/s = 200 GeV (right) 
calculated with Pythia. Pions, kaons and protons/neutrons 
are seperated by line colour, while their charge state [negative 
(long dashed), neutral (solid), positive (dashed)] is indicated 
by the line style. Neutral anti-particles are displayed with a 
thin line. The sum of all charged particles is depicted by the 
thick dotted curve. Note that the spectra of all pions are sim- 
ilar and the spectra of neutral (anti) kaons and nucleons are 
similar to that of their positive (negative) charged partners. 



lisions at these energies for midrapidity are given by 
respecting the correct isospin average corresponding to 
0.16 • pp + 0.36 • nn + 0.48 • pn. It is worth noting, that 
these distributions do not follow a true power law be- 
haviour, as can be seen in Fig. |3 where at p± ~ 2 GeV 
and ~ 6 GeV significant changes in the slope occur. 

Some special care has to be taken (for example mix- 
ing of calculations with non-weighted and some with 
reweighted events) to expand the Pythia Monte-Carlo 
calculations up to values of p± larger than 10 GeV with 
best statistics. Naive Monte-Carlo implies, that at p± w 
6 GeV reasonable computation time ends and that one 
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would have to increase the number of events by approxi- 
mately a factor of 10 every time one wants to go further 
in the p± spectrum by an other step of 2 GeV, as one can 
see in Fig. |3J 

B. Folding Matrices 

The collision of a high p± particle with particles of the 
medium is implemented according eq. Q via a folding 
matrix gij(p± ,px)- Since a unique description of scat- 
terings at these low energies (cf. Fig. is not available, 
model descriptions have to be choosen. 

In our first description A we calculate the folding ma- 
trices with Monte-Carlo calculations using the package 
Fritiof [2(j to simulate collisions with a p-meson. Since 
it is known, that within Fritiof the relative weight- 
ing between elastic and inelastic scatterings is not in 
full agreement with data, we therefore choose a constant 
weighting according the ratio 15:85. (This is given by 
the averaging of p + p-, it + tt-, K + ir cross sections 
in a pomeron/reggeon exchange picture for y/s > 2 GeV 
and also following simple estimates in a constituent quark 
model.) 

As a second description B we use purely elastic scat- 
tering on a pion, where the angular distribution of the 
outgoing particles in the cm system is assumed to be 
isotropic, which seems to be a at a first glance a rather 
strong assumption. 

In the description A of the collisions we are forced 
to extrapolate the Fritiof prescription downwards to 
energies, where it is strictly not fully applicable any- 
more (for values below an initial transverse momen- 
tum of p± ~ 3 GeV we are in the resonance region, 
cf. Fig. 2J. Therefore (or nevertheless) we stop the cal- 
culations at a minimal initial transverse momentum of 
about p"" n = 1 GeV. Boundary problems applying the 
folding procedure acccording eq. Q due to this cut could 
be easily cured, if one implements a more detailed weight- 
ing of inelastic and elastic cross sections, which could al- 
ready allow to bend the calculations downwards to zero 
transverse momentum, or if one would be able to intro- 
duce a consistent picture of scattering of hadrons from 
very low (resonance region) to very high momenta (hard 
QCD processes). 

The two different model descriptions show some sig- 
nificant differences, if one looks at the dependence of the 
final transverse momentum p± for a given initial trans- 
verse momentum p± . But not only the different models, 
but also different final particles behave differently. This 
is shown in Fig. where as an example the convolution 

matrizes g^+/ p + and gf^+i v + are shown for the model 
descritpions A and B for the fixed value of the initial 
transverse momentum p± = 5 GeV. 

While all the elastic distributions (Fig. 0] (bottom)) do 
follow a flat behaviour, where the threshold values are 
dictated by p± and the corresponding masses of the par- 
ticles, the inelastic distributions (Fig. 0] (top)) do possess 
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FIG. 4: The folding matrices g according (in)elastic scatter- 
ing on a p (top) or elastic scattering on a it (bottom) for a 
7r + (left) and a p (right) as outgoing particle for a fixed value 
of the initial transverse momentum p± — 5 GeV as function 
of the resulting transverse momentum p± . The line styles are 
choosen as in Fig. |3J indicating the different incoming parti- 
cles. 'Negative' values for transverse momenta p± stand for 
scattering in the opposite direction with positive momenta. 




FIG. 5: As Fig. 0] but now as function of the initial trans- 
verse momentum p± for a given resulting transverse momen- 
tum p± — 5 GeV. The connected dots stand for the slope of 
/(ft++h-)/2(P±o) at Vs = 200 GeV. 



a much richer structure. First they show, as expected, 
a strong enhancement at low p± values. In addition, 
these curves also depict a double peak structure for the 
largest p± -values in the forward scattering region. The 
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lower peak at p± — pj_ a is the simple one due to elastic 
scattering, while the second and rather astonishing peak 
for p± > p± is due to a large partial cross section for 
X + p — > X + 7r implemented in Fritiof (X denotes 
here either the tt + in the upper left figure or the proton 
in the upper right one of Fig. 0J). This peak is located 
at (p x - p± ) ^ (m 2 p - ml)(2m p ) w m p /2. Physically, 
for such collision events the leading particles do become 
slightly accelerated. One might argue, that the likelihood 
of such contributions are completely overestimated, so 
that Fritiof should also fail in this part of the descrip- 
tion. On the other hand, by including them, their net 
effect in our later discussion would be that one actually 
underestimates the emerging depletion of the transverse 
momentum spectrum. We have decided to omit these 
contributions from the inelastic part by hand compared 
to the basic Fritiof concept. 

Considering eq. (0J, one realizes the fact that the be- 
haviour of the folding matrices for a fixed value of the 
resulting transverse momentum p± as a function of the 
initial transverse momentum p± is of crucial importance. 
The corresponding figures for the examples from above 
are depicted in Fig. [5] 

While the curves for the elastic scattering on a pion 
are smoothly falling for p± > p± (or the corresponding 
threshold value), the curves for the scattering on a p are 
nearly constant, but also showing the double peak struc- 
ture at p± = p± and p± > p± , as mentioned above. 
For larger p± the behaviour of all these folding curves 
gets less and less important for the final spectrum, since 
one has to fold them with the ininital distribution, whose 
steep slope is also indicated in Fig. [5] This means that 
the correct way of looking at the stopping is the function 
gij(j>± = const) instead of gij(px = const). Looking 
at the potential size of the energy loss (at a partonic 
or a hadronic level), the value p± —p± of gij(p± ,P±) is 
strongly skewed toward small values by the steeply falling 
distribution /°(p± ). Only the least stopped particles 
(i.e. particles with pj_ only somewhat larger than p±) in 
gij contribute to the folding integral eq. (@J. Simple esti- 
mates give, that values of with p± < \px a contribute 
only at a one percent level. Although appealing from the 
various shapes of the folding matrix g^j , assumptions like 
AE = Ap± = (p± — p± ) ~ const ■ p± for an averaged 
energy loss are much too oversimplistic. 

These findings are in complete agreement with the con- 
siderations of |7j , where the energy loss of (very energetic) 
partons in partonic matter due to gluon radiation, calcu- 
lated via perturbative QCD, is expressed via a function 
Z?(e), which has to be integrated in a sensible way with 
e as a measure of the energy shift. 

Inspecting further Fig. El the p± > p± peak of the 
Fritiof results will lead to some (maybe unphysical) 
"unstopping" , as already emphasized. As those contri- 
butions are indeed then weighted with the much larger 
distribution f(p± ) at p± < p±, they would indeed be 
of more crucial importance than naively expected. In- 
cluding those contributions would roughly weaken the 



deacceleration effect in the spectrum by a factor of 2 for 
Pj_ > 3 GeV and lead to a zero net effect at p± ~ 2 GeV. 
As already discussed, we will exclude these contributions 
with px < p±- 



III. QUENCHING BY MULTIPLE COLLISIONS 
AND COMPARISON WITH DATA 

Employing the appropriate scaling of binary collisions 
(N coU = 945, 905, 19, 3.7 for the centralities 0. . . 5%, 
0...10%, 60... 80% and 80... 92% respectively), the 
Pythia calculations are in perfect agreement with the 
(at that time preliminary) PHENIX data 25] for charged 
hadrons and for neutral pions in the case of peripheral 
collisions, as can be seen in Figs and We can thus 
conclude that our Pythia calculation with the correct 
extrapolation provides an accurate basis for the experi- 
mental comparison. 

Looking at the most central data, one recognizes the 
celebrated significant discrepancy between our scaled 
Pythia results and the data. The data resides below the 
theoretical curve for p± < 1.5 GeV and for p± > 2.5 GeV. 
Only the latter case for p± > 2 GeV we now want to 
adress by means of our scenario of possible (multiple) 
collisions of the jet-like (pre-)hadrons with the 'soft' 
hadrons sitting in the bulk matter of the fireball. There 
is some general agreement that the hadron spectra below 
p± ~ 2 GeV should be dominated by soft physics and 
hydrodynamical expansion with some appreciable collec- 
tive outward flow. In addition, in this lower momentum 
regime our basic assumption of a hadronic target being 
approximately at rest becomes no more valid. Here not 
only energy loss, but probably also some collisional en- 
ergy gain (on a partonic and/or hadronic level) and also 
hadronic particle annihilation and (re-) absorption occurs 
which influences the shape of the spectrum. 

Refering to eq. I0J and inspecting the two Figs. 0] and 
[SI one recognizes that the depletion in the spectrum by 
one collision should be stronger for description A with 
(in)elastic collision on a p as the soft target particle, as 
one would also intuitively expect. This is so as the fold- 
ings matrices g^ for pj_ > p± are somewhat smaller 
for this first collision picture, though there is the same 
strength at p± a ~ p± . In Figs and [7| the three solid 
lines depicts the spectrum (for p± > 2 GeV) for the most 
central region, when the particles have suffered on av- 
erage 1,2,3 (in)elastic hadronic collisions on a p. The 
dotted lines in these figures indicate the effect for elas- 
tic scattering on a i. Indeed, for the inelastic case the 
overall depletion is somewhat stronger. Refering to our 
earlier remarks, we only do show the modifications for 
p± > 2 GeV. One clearly recognizes the potential im- 
portance of these possible final state interactions. For 
p± > 2 GeV an average loss with (L/X) = 1 ■ • • 2 colli- 
sions are appropriate to explain the present experimental 
data both for charged hadrons and neutral pions. 

We remark, that the potential energy loss by final state 
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FIG. 6: Comparison of the calculations at midrapidity with 
the (at that time preliminary) PHENIX data |25| (but see 
also (J) on charged hadrons for three different centralities. 
The thin dashed lines indicate the results from Pythia scaled 
with a corresponding binary collision number N co ii- The solid 
lines depicts the spectra for the most central region, where the 
particles have suffered on average (L/X) = 1, 2, 3 (in)elastic 
hadronic collisions on a p, while the dotted lines indicate the 
same effect for elastic scattering on a ir. These modifications 
are only shown for p± > 2 GeV. 
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FIG. 7: As Fig. El but here a comparison of the present 
calculations with the PHENIX data [23 ] on 7ro is depicted. 



hadronic reactions is in the same range as considered in 
the various discussions for the jets in possible deconfincd 
matter. 

In Figs. and 03 we now show calculations for vari- 
ous number of collisions at a ^/snn = 200 GeV up to 
p± = 10 GeV. Our results might serve also as a guide- 
line and prediction concerning the present analysis of the 
200 GeV runs at RHIC 3]. In principle, the results are 
completely similar to the one carried out at the lower 
RHIC energies. Fig. [S] shows the resulting p± spectra af- 
ter 1,2,3 collisions according to the two collision pictures. 
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FIG. 8: Resulting px-spectra of charged hadrons at midra- 
pidity for ^/s = 200 GeV (solid black line) and for (L/X) = 
1, 2, 3 (top to bottom) collisions according (in)elastic scatter- 
ing on a p (blue) or elastic scattering on a ir (green- light). 
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FIG. 9: The suppression factor R(p±) of charged hadrons 
at midrapidity for yfs = 200 GeV for (L/X) = 1, 2, 3 (top 
to bottom) collisions according (in)elastic scattering on a p 
(blue) or elastic scattering on a 7r (green). 



As expected, the modification of the spectrum is similar 
to the one discussed for the situation at lower energies 
of y/s = 130 AGeV. Thus we can conclude that hadronic 
final state collisions are efficient in degrading the trans- 
verse energy of jet-like hadronic particles. 

Fig- El gives the resulting suppression factor R(p±) = 
= J modificd for collisions from Fig. EI Here the effect 

litial 

With 



dN/dp ± \ 

of the different collision schemes is more visible 



one to two inelastic collisions one will have a suppression 
of 50 to 80 percent for momenta p± > 4 GeV. 



IV. SUMMARY 

We have motivated that most of the (pre-)hadrons 
stemming from a jet should still materialize in the dense 
system for transverse momenta up to 10 GeV. The late 
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hadronic final state interactions with the bulk of co- 
movers have a clear and nonvanishing effect in suppress- 
ing the spectrum. This is so for (in)elastic reactions on 
a p, which we have modelled via the FRITIOF scheme, 
and also for (isotropic) elastic scattering on a n. On the 
average one up to two such interactions should already 
be enough to explain quantitatively the RHIC results. 

Although this comparison of the depletion of the spec- 
tra with real data is intriguing, one should be aware that 
our calculations are at best semi-quantitative. This is so 
because of a couple of simplifying assumptions invoked 
in the present study. We have neglected all possible de- 
pendencies of the cross sections of the colliding "hard" 
and "soft" (target) particles, but do assume that a fixed 
number of collisions has taken place, following a simple- 
minded opacity expansion. In particular, possible res- 
onance contribution for the elastic scattering could en- 
hance the effect of momentum degradation considerably. 
Second we made the further assumption, that the bulk 
hadronic targets are taken at rest. This is certainly not 
the case especially if one wants to quantitatively discuss 
the modifications in the spectra around p± = 1 . . . 3 GeV, 
where both collective and 'thermal' dynamics should play 
a considerable role. Energy gain and loss mechanisms 
should both be at work and are probably competing in 
this momentum range. On a partonic level such a picture 
has been discussed in [26| , but this should be rather obvi- 
ous if one thinks of a kinetic cascade-type description ei- 
ther on a partonic or hadronic description of the dynam- 
ics happening inside the fireball. Hence, the marriage of 
a transport description for the 'soft' matter by some stan- 
dard ultrarelativistic heavy ion cascade together with the 
present ideas carried out in this study is one possible way 
to further proceed in order to obtain more quantitative 
results. We will leave this as a future option. A first 
and rather exploratory investigation in this direction has 
already been mentioned some time ago |27| . 

Our emphasis is, however, at this stage to point out 
that, in order to draw some deductions for possible QCD 
effects of a deconfined QGP phase on the materializing 
jets, one has to disentangle these from the here investi- 
gated final state interactions, before definite conclusions 
on the importance of a potential dense partonic phase, or 
any other effect, can really quantitatively be drawn. This 
is certainly true for the case of the jets to be observed at 
RHIC. (At future LHC experiments observable jets with 



much higher p± could actually materialize outside the 
fireball.) 

Combining specific energy loss fomulae with fragmen- 
tation functions of quarks and gluons into hadrons can 
therefore only be understood as an extremely "effective" 
approach, since no conclusive answer about the relative 
strength of hadronic final state interactions can be ex- 
tracted. The possible disentanglement, though, is very 
challenging and delicate. One possible way could be 
the precise measurements of p± -distributions of charmed 
particles like the .D-mesons [28|. According to a recent 
investigation [2^| heavy quarks should loose consider- 
able less energy than light quarks or gluons (although 
an earlier study comes to other conclusions [30f ). On the 
other hand the (elastic and inelastic) D-meson hadronic 
cross section should not be much smaller than that of 
light quark mesons. Hence, if the .D-mesons should 
not show any significant p± -suppression at moderate 
p± ~ 3 ... 10 GeV, this would then favour the partonic 
'pQCD' jet quenching scenario. At present, the PHENIX 
data [31j up to p± values of 3 GeV still suffer from too 
large statistical errors and are compatible with full or no 
quenching |3^ |. In any case, the D-meson spectra will 
be an important experimental challenge for future RHIC 
undertakings. Some other valuable information on the 
right "initial distribution" will be learned by d + Au ex- 
periments |12| ], eg. the size of possible shadowing and the 
Cronin effect. However, like the medium-induced gluon 
radiation the here discussed final state interactions with 
the late stage hadronic matter is a pure heavy ion ef- 
fect. Here the p^-dependence of the 'suppression' pat- 
tern might reveal possible differences, as for higher p± 
the time of formation for the hadronic wavefunction of 
the leading hadrons of a jet and the colour transparency 
effect should suppress these final state interactions. But 
this is speculation at present. 
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